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Abstract. Most models for the central engine of gamma-ray bursts involve a stellar mass black hole surrounded 
by a thick disk formed after the merging of a system of compact objects or the collapse of a massive star. Energy 
released from the accretion of disk material by the black hole or from the rotation of the hole itself extracted 
by the Blandford-Znajek mechanism powers a relativistic wind along the system axis. Lorentz factors of several 
hundreds are needed to solve the compactness problem in the wind which implies the injection of a tremendous 
power into a very small amount of matter. The Blandford-Znajek mechanism, where the outflow follows magnetic 
field lines anchored to the black hole is probably the best way to prevent baryonic pollution and can even initially 
produce a purely leptonic wind. In this paper we rather study the wind emitted from the inner part of the disk 
where the risk of baryonic pollution is much larger since the outflow originates from high density regions. We show 
that the baryonic load of this wind sensitively depends on the disk temperature and magnetic field geometry and 
that the outflow can become ultra-relativistic (Lorentz factor T > 100) under quite restrictive conditions only. 
Conversely, if F remains of the order of unity the dense wind emitted from the inner disk could help to confine 
the central jet but may also represent a source of baryon contamination for the Blandford-Znajek mechanism. 



Key words. Gamma rays: bursts 
Relativity 



Accretion: accretion disks - Magnetohydrodynamics (MHD) - Neutrinos 



1. Introduction 



The di znwry of the first optical rxnrrrfcerprixts Lu gaiiiina- 
ray bursts (hereafter GRBs) hi 1997 (van Paradijs el al 



et al. 1992; Meszaros 



systems) ( 


Narayan 


Mochkovitch et al. 



1993 ) or massive stars which collapse to a black hole 



1997) has shown that most (if not all) GRBs are located 
at cosmological distances. About 20 redshifts have now 
been measured from z — 0.43 for GRB 990712 to z = 4.5 
for GRB 000131 (with however the peculiar case of 
GRB 980425 which appears to be associated to a near by 
type Ic supernova at z = 0.01; Galama et al. ( 1998j ) ) . 
The energy radiated by these cosmological GRBs in the 
BATSIf} 



-fcrr 



range (20 - 1000 keV) goes hum 510^ eigs 
GRB 9 70228 and GRB 980013 (at z - 0.095 and 1.090) 

to 2 10 54 ergs for GRB 990123 (at z = 1.6) assuming 
isotropic emission. After correction for beaming, the true 
energy output appears to be less scattered, clustered 



(collapsars) ( Woosley 1993 ; Paczynski 1998 ). In all these 
cases, the resulting configuration is a stellar mass black 
hole surrounded by a thick torus made of stellar debris 
or of infalling stellar material partially supported by 
centrifugal forces. The location of the detected optical 
counterparts well inside their host galaxies and often 
associated with regions of star formation appears to favor 
the collapsar scenario ( Paczynski 1998| ; Owens et al 



around E~ ~ 5 10 ergs (Djorgovski et al. 2001). Among 



the sources which have been proposed to explain such 
a huge release of energy in a time of seconds, the most 
popular are mergers of compact objects (neutron star 
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1998; Klose et al. 200C). Double neutron star or neutron 
star - black hole mergers should generally be observed at 
the periphery of the host galaxy due to the long delay 
before coalescence and the large velocity imparted to 
these systems by two successive supernova explosions. 
They can however still be invoked in the case of shorts 
bursts, for which no optical counterpart has been detected. 

If black hole + thick disk configurations are indeed at 
the origin of GRBs, the released energy will ultimately 
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come from the accretion of disk material by the black hole 
or from the rotational energy of the hole itself extracted 



by the Blandford-Znajek mechanism (Blandford & Znajek 
19770~|ln a first step, the energy must be injected into 
a relativistic wind whose existence has been directly 
inferred from the observations of radio scintillation in 



GRB 970508 ( |Waxman et al. 1998| ) and which is also 
needed to avoid photon-photon annihilation and the 



resulting compactness problem (Baring & Harding 1997). 
The second step consists of the conversion of a fraction 
of the wind kinetic energy into gamma-rays vi a the 



format ion of shoc ks, probably inside the wind itself ( |Rces 



& Mes ;aros 1994 ). These internal shocks can be expected 
if the source generate a highly non uniform distribution 
of Lorentz factor so that rapid layers of the wind will 
catch up with slower ones at large relative velocities. In 
the last step, the wind is decelerated when it interacts 
with the interstellar or circumstellar medium and the 
resulting (external) shock is responsible for the afterglow 



observed in the X-ray, optical and radio bands (Meszaros 
& Ree j 19971) . 



The physics of the afterglow is probably the best 
understood since most afterglow properties can be 
interpreted in terms of solutions of the relativistic Sedov 
problem (Blandford & McKcc 1976) with synchrotron 
emission behind the shock (Sari et al. 199S). Models 



taking into account the beam geometry of the flow 
( [Rhoads 1997 ) or different kinds of burst environments 



(constant density medium or stellar wind; Chevalier & Li 
(2000|)] can be constructed to explain observed breaks in 
the lightcurves or the evolution of some radio afterglows. 

More problems remain concerning the generation of 
gamma-rays in the relativistic wind. Instead of internal 
shocks, gamma-rays can also be emitted during the early 
evolution of the external shock with however the difficulty 
to explain in this case the highly variable temporal 
profiles of observed burst s ( Sari et al. (199§| ) see however 
Dernier fc Mitman (1999 )). Models of bursts produced by 
internal shocks are found to be in reasonable agreement 
with t he observation s ( Kobayashi et al. 1997 ; Daigne 



& Mo chkovitch 199S , 2000 ) even if the process which 
produces the gamma-rays - synchrotron radiation or/and 
comptonization - remains uncertain due to problems 
encountered in fitting the low energy part of the spectra 
( |Preccc ct al. 1998j ). 



The origin of the relativistic wind is the more complex 
of the three steps. Several proposals have been made to ex- 
plain the generation of this wind but few detailed calcula- 
tions have been performed. If the burst energy comes from 
matter accretion by the black hole, the annihilation of 
neutrino-antineutrino pairs emitted by the hot disk could 
be a way to inject energy along the system axis, in a re- 
gion which can be expected to be depleted in baryons due 
to the effect of centrifuga l forces ( Meszaros fc Rees 1992a ; 
Mochkovitch ct al. 1993| ). The low efficiency of this pro- 



cess requires high neutrino luminosities and therefore high 
accretion rates. In the merger case this can be achieved for 
short accretion timescales (Ruffcrt & Janka 1999) and may 
explain short bursts. Conversely in the collapsar scenario, 
the larger mass reservoir allows the system to maintain 
high accretion rates for a longer time and could then also 



produce long bursts (MacFadyen & Woosley 1999). 

Another possibility is to suppose that disk energy is ex- 
tracted by a magnetic field amplified by differential rota- 
tion up to very large values (B <; 10 15 G). A magnetically 
driven wind could then be emitted from the disk with a 
fraction of the Poynting flux being e ventually transferred 
to matter (|Blandford fc Payne 19821) . Such a mechanism 
operates in many classes of astrophysical objects (from T 
Tauri stars to AGN) but it is far from clear that it can 
work in the context of GRBs where final Lorentz factors 
of several hundreds are required. In a different version of 
the same idea, an early conversion of magnetic into ther- 
mal energy could occur through the reconnection of field 
lines above the plane of the disk in a region of rather low 
density ( Narayan et al. 1992 ). 

An alternative to accretion is to directly extract the 
rotational energy of the black hole via the Blandford- 
Znajek mechanism. The available power then depends on 
the rotation parameter a of the hole and on the intensity 
of the magnetic field pervading the horizon. If B <; 10 15 G 
and a ~ 1, the power available from the Blandford-Znajek 
mechanism can be larger than 10 52 erg.s^ 1 with a very 
limited contamination by baryons at the source since the 
field lines which guide the outflow are anchored to the 
black hole. 

In this paper we rather concentrate on the wind which 
is emitted from the inner disk. We want to identify the 
key parameters which control its baryonic load and check 
whether it can reach large Lorentz factors or remains non 
relativistic. Our approach will be oversimplified in com- 
parison to the complexity of the real problem so that our 
conclusions have to be considered as indicative only. In 
Sect. H we briefly discuss the structure of the disk + black 
hole configurations which are produced by NS + NS or 
NS + BH mergers and in the collapsar scenario. We write 
in Sect. || the equations which govern the wind dynamics 
from the disk to the sonic point. They are solved in Sect. |^ 
to obtain the mass loss rate and the terminal Lorentz fac- 
tor of the wind is estimated in Sect. |^. Our results are 
discussed in Sect. || which is also the conclusion. 

2. The structure of the disk 

Before writing the wind equations in the next section we 
first describe the black hole + disk configurations which 
are obtained in the case of the three most discussed GRB 
sources: NS + NS or BH + NS mergers and collapsars. 
We want to compare the black hole and disk masses, the 
hole rotation parameter a, estimate which fraction of the 
disk is optically thick to neutrinos and obtain the disk 
temperature. 
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Fig. 1. Equilibrium models for thick accreting disks orbi- 
ting stellar mass black holes. The disks are supported by 
the pressure of relativistic degenerate electrons and the 
specific angular momentum is increasing outwards as r 2 ; 
case 1) Schwarzschild black hole of mass Mbh = 2.5 Mq 
and disk mass Md = 0.3 Mq; 2) same Mbh and Md 
and Kerr black hole with a = 0.8; 3) Mbh = 5 Mq, 
Mb = 0.5 Mq and a = 0.4. In each model thick lines 
in the equatorial plane indicate the extension of the zone 
which is optically thick to neutrinos for T v = 2 (left) and 
3 (right) MeV. 



2.1. Mergers : 

The coalescence of neutron stars has been studied 
by several groups mostly in the Newtonian and post- 



Newtonian a 


□proximation (Davies et al. 1994; 


Ruffcrt 


3t al. 1996 




Calder et al. 1999 




Rosswog et al. 2000; 


Faber et al. 200C) The resulting 


merged object obtained 



from two neutron stars of 1.4 Mq is made of a dense 
central core of ~ 2.5 Mq in quasi-uniform rotation 
surrounded by a differentially rotating disk of ~ 0.3 Mq. 
In relativistic calculations, a black hole is directly formed 
during the merging event if the total mass of the system 
is 30 to 70% larger than the maximum rest mass of an 



isolated neutron star (Oohara & Nakamura 1997). Due to 
the large angular momentum of the merger the rotation 
parameter a = Jc/GM^ K of the newly formed black hole 
can easily exceed 0.5. 

In the case of BH + NS mergers, the masses of the 
disk and black hole are in average larger than in NS + 
NS mergers. The hole rotation parameter depends on 
the fraction of neutron star materi al accreted by the 
black hole during the merging event. Janka et al. (1999 ) 
obtain disk masses between 0.3 and 0.7 Mq and rotation 
parameter between 0.1 and 0.5 for different assumptions 
regarding the relative masses and spins of the neutron 
star and the black hole. 

To estimate the transparency of the disks to neutri- 
nos we have computed the structure of the merged object 
with the self-consistent field method originally developed 
by Ostriker & Mark (1968) using the approach of Hachisv, 
|(1986| ) to solve the Poisson equation. The equation of state 
in the disk corresponds to an ideal gas of ultra-relativistic 
electrons and the distribution of specific angular momen- 
tum is j(r) oc r 2 in reasonable agreement with the results 



of numerical simulations (Ruffcrt & Janka 1999). We con- 
struct accreting disks i.e. disks with a cusp in the equa- 
torial plane where the total gravitational + centrifugal 
force is zero. If a = 0, we model the black hole with the 
Paczynski-Wiita potential while for a Kerr black hole we 
use the Novikov potential ( Artemova et al. 199G| ) 



$bhM = - 



GM 



BH 



(/3 " l)r fc 



0-1 



with 



?7l 



(l + Vl-a 2 ) 



and 



(3=^1 



(1) 



(2) 



where r g = GMbh/c is the gravitational radius, is 
the horizon radius and r ms the radius of the last stable 
orbit. We have considered three cases: Mbh = 2.5 Mq, 
M d = 0.3 M Q with a = (case 1) and a = 0.8 (case 2); 
M BH = 5 Mq, M d = 0.5 Mq with a = 0.4 (case 3). Cases 
1 and 2 are representative of NS+NS mergers while case 3 
corresponds to a BH+NS merger. The obtained disk struc- 
tures are represented in Fig. [y. Compared to case 1, the 
increase of a in case 2 leads to a denser (/9 max = 4.2 10 12 
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instead of 3.3 10 11 g.cm -3 ) and more compact disk (the 
radius at maximum density r max and the external radius 
r cxt are 58 and 279 km in case 1, 25 and 117 km in case 2). 
In case 3, the maximum density is /0 max = 4.3 10 11 g.cm~ 3 
at fmax = 80.4 km and the external radius is r cxt = 261 
km. We have computed the optical thickness of these disks 
to electron neutrinos in the vertical direction 

r+oo 

, ,. = / n Ve p dz (3) 



where 



3.8 10 



-19 



cm .g 



(4) 



1 MeV 

is the neutrino opacity ( Duncan et al. 1986| ). The radial 
extension of the optically thick region has been repre- 
sented in Fig. |l| for two temperatures, T v = 2 and 3 MeV. 
The optical thickness becomes larger than unity at 3 - 5 
r/i and the disk remains opaque out to r out ^ 10 rn- 

When the disk is optically thick we use the following 
expression for the neutrinosphere temperature 

1/4 

TAn = T{-\ I 1 - ^% \ (5) 




where rj n is the disk internal radius and T* is the 
temperature at a reference radius (we take below 
= 4 r/j). This kind of behavior is expected in the case 
of a geometrically thin disk but is certainly a very rough 
approximation in the case of a thick disk. 



2.2. Collapsars : 

In collapsars the accretion flow toward the black hole has 
been studied in detail by Popham et al. (1999 ). The in- 
ner disk is fed by material from the collapsing stellar en- 
veloppe. It is less dense than in the merger c ase (typical 
densities are 10 8 g/cm 3 ) and is optically thin (MacFadycn 
fc Woosley 19991) . We follow |Popham et al. (19991 ) to ob- 
tain an analytical expression for the temperature from the 
balance between the dissipated energy and neutrino losses 

QT e = -v Vl\ (6) 



where q e N is the cooling rate per unit mass due to the 
emission of neutrinos by nucleons (Q = 1.4 10 18 with T in 
MeV); q e N is the dominant cooling contribution as long 
as the disk temperature does not exceed about 10 MeV. 
We adopt an a-prescription for the disk viscosity 

v = aHvs (7) 

where H — is the disk half thickness and v 8 is the 
sound velocity. If the perfect gas contribution dominates 
in the disk then, v 2 ~ (p being the average molecular 
weight) and the disk temperature is given by 

-0.3 





Fig. 2. Schematic view of the system geometry. The disk 
is assumed to be geometrically thin and the magnetic field 
lines make an angle 8 with the disk. The outflow follow 
the field lines and the position of a fluid element along a 
line anchored at a radial distance r from the black hole is 
represented by y = i/r. 

3. Dynamics of the wind from the disk to the 
sonic point 

3.1. Wind equations 

To study the dynamics of the wind we make a number 
of simplifying assumptions. We first suppose that the 
inner disk is geometrically thin which is certainly wrong 
in the merger case and probably a poor approximation 
in the collapsar scenario. Wind material is guided along 
magnetic field lines for which we adopt the simplest 
possible geometry: close to the disk the field is poloidal, 
made of straight lines making an angle 9(r) with the 
plane of the disk (r being the distance from the foot 
of the line to the disk axis; see Fig. |^). Since we limit 
our study to the part of the flow below the sonic point 
(which is sufficient to obtain the mass loss rate) we use 
non relativistic equations (v/c <C 1) but we adopt a 
Paczynski-Wiita potential for the black hole. We also 
assume that the wind has reached a stationary regime. 
Clearly, a realistic description would imply a thick disc, 
a complicated field geometry and a time-dependent wind 
dynamics but we believe that the toy model presented 
in this paper remains able to identify the main physical 
processes which affects the baryonic pollution of the wind. 

We write the three flow equations in a frame corotating 
with the foot of the line: 

— Conservation of mass : 
pvs(y) = m , 

— Eulcr equation 

i(y)r - 



dv 



1 dP 

dy p dy 

Energy equation : 



(8) 



de 
dy 



q(y)r 



P_dp 

P 2 dy 



(9) 
(10) 

(11) 
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where y — £/r, I being the distance along the field line 
(£ = in the plane of the disk); e is the specific internal 
energy, 7(2/) and q(y) are respectively the total accelera- 
tion (gravitational + centrifugal) and the power deposited 
per unit mass in wind material due to neutrino heating 
(and cooling), viscous or ohmic dissipation and magnetic 
reconnection. Because the field and stream lines are coin- 
cident the function s(y) can be simply related to the field 
geometry. We get 



s (y) = l + ay + by 2 



with 

a = cos 9 - 
and 

b = — cos ( 



(10 



smf 



suit 



dlogr 



(16 



(12) 



(13) 



(14) 



dlogr 

Finally, rh is the mass loss rate per unit surface of the 
disk. The acceleration j(y) is derived from the potential 

dy 



i{y)r = 

where 
$BH(y) 



$(y) = $BH(y) + $c(y) 



(15) 



GM- 



BH 



r^Jy 1 + 2ycos9 + 1 - r h 



1 



2y cos 9 + 1 - 1 



(16) 



is the black hole potential (x — r/r^) and 
$ c (y) = -\tt 2 {r)r 2 (l + ycos9) 2 



c 

T 



(x-iy 



r(l + ycos9) 2 



(17) 



is the centrifugal potential. There is a critical angle 9 cr 
below which the acceleration is always positive so that 
matter (even at zero temperature) can escape freely from 
the disk without being confined in a potential well. In new- 
tonian gravity 6 cr = 60° ( Blandford fc Payne 1982 ) but 
9 CI slowly decreases when the Paczynski-Wiita potential 
is used, from about 63° at x — 3 to 60° at large radial 
distances. 



3.2. Equation of state 

Our equation of state includes nucleons, rela tivistic elec- 
trons and p ositrons, and photons. Following Bcthe et al 



(198C ) and Bcthe (1993 ) the contribution of relativistic 



particles is given by 



Pr = 



(fcT) 4 /llTT 



12(?ic) 3 V 15 
(1.26 + 0.3577 2 



0.017r/ 4 )10 26 T^y dyne.cm -2 



(18) 



where r\ = p, e /kT, \x e being the electron chemical poten- 
tial. Nucleons behave as an ideal gas of pressure 



P N = -?—kT 

m N 



(19) 



where mjy is the nucleon mass and the density is obtained 
from 



3 \hc J 



7.2 10 7 T^Yri (77+^ 1 -..en, 



-3 



(20) 



The number of electrons per nucleon Y e should be com- 
puted from the rates of neutrino capture and emission by 
nucleons. In practice, we do not perform this calculation 
and simply adopt a constant Y e = 0.5. Finally, the specific 
internal energy is 



(3P r + 3/2P N ) 



(21) 



so that the three thermodynamic quantities P, p and e 
that appear in the flow equations can be expressed in 
terms of T and 77. 

3.3. The sonic point 

To solve the flow Eqs (|9|-|ll|) we first derive Eq. (||) and 
express the thermodynamic variables p, P and e as func- 
tions of T and 77 to obtain a linear system for log v, log T 
and log 77 



dlogv , g rflogT | dlogr? 



dy 
2 d log v 



dy 



P 



dy 

^PdlogT 

P dy 
dlogT 



dy 

^ log 77 



dy 



dy 

C(v)e - A{rf) 



dlogsjy) 
dy 

y(y] 

d log 77 



P 



dy 



<i(y)r , 

(22) 



where we have used the equation of state to get 



<91ogp 



dlogT 
<91ogP 



<91ogp 



dlogT 
dloge 



d log 77 
dlogP 



aiog-r 



1 



d log 77 
dloge 



d log 77 



--B( V ) , 



(23) 



The three derivatives of v, T and 77 can then be written as 

dlogv = Fi(t/,t;,T, 77) 

dy A 

dlogT = F 2 (y : v,T,7 1 ) 

dy A 

dlogrj _ F 3 (y,v,T,r}) 

dy A 



(24) 
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where A(v,T,r]) is the determinant of the system. It is 
equal to zero at the sonic point which gives 



v 2 = 



(B-4<7)-(7-l)(3fl-4A) 



(A - 3C) 



P 
P 



(25) 



In the case of an optically thick disk we make the extreme 
assumption that all the heating which is not due to neu- 
trino processes takes place below the neutrinosphcre. The 
power q injected into the wind is then restricted to the 
neutrino contributions 



where 7 — 1 = Pipe. If the equation of state is dominated 
by the contribution of relativistic (resp. non relativis- 
tic) particles 7 — 1 = 1/3 (resp. 2/3) and?; 2 = |(resp. §)^. 

Since the three functions v, T and 77 must remain re- 
gular everywhere in the wind the numerators in Eq. ( p4[ ) 
must be zero at the sonic point 

Fi(y 8> v s , T s , T] s )=F 2 (y s , v s ,T s , T] s ) = F 3 (y s ,v s , T s , r) s )=0 

(26) 

which yields a unique relation among the parameters at 
the sonic point 



Qv — Qi/N + <lve + ivv — q e N — ie+e 



(30) 

where the different terms in Eq. (^0|) respectively cor- 
respond to neutrino capture on free nucleons, scatter- 
ing on electrons and positrons, neutrino- antineutrino- 
annihilation (heating) neutrino emission by nucleons and 
annihilation of electron-positron pairs (cooling). With our 
assumption that Y e = Y p = Y n = 0.5 the electron neu- 
trino and antincutrino temperatures are identical and the 
heating by capture on free nucleons takes the simple form 
( |Qian fc Woosley 1996Q 

L) , T (4.Y 



Qi/N 



dlogs(y) 



dy 



q{y a )r ( P_ AA-3B 
pe A-3C 



disk 



dft 



= 



1.4 10 18 / T, 



(27) 



disk 



dn 

47T 



erg.g 



(31) 



the parenthesis being equal to 1/3 (resp. 2/3) for relativis- 
tic (resp. non relativistic) particles. Since (?) the last term 
in Eq. ( pTj ) is smaller than the two others in most cases 
of interest and {it} the derivative of s{y) is positive (if the 
inclination angle 9{r) of the field lines decreases with in- 
creasing distance to the axis) it can be seen that the sonic 
point is located in the region where 7(1/) is negative, i.e. 
below yi where 7(2/1) = 0. It appears in practice that y s 



7 <rT* ■ 

— '- is the 

16 7T 

(e™ ) are the 



where the temperature is in MeV, I Ve = Ip e 
neutrino (antineutrino) intensity and (e™ ) 
nth energy moments of the neutrino Fermi distribution. 
The integral is performed over the disk surface and dft is 
the solid angle of a surface element as seen from a point 
of coordinate y on the field line. All relativistic effects 
on the neutrinos (bending of trajectories, gravitational 
and Doppler shifts) have been neglected even if they can 
lead to appreciable corrections in the results (Jaroszynski 



Vi 



3.4. Heating and cooling sources 

Several sources can contribute to the injection of energy 
in wind material: viscous or ohmic dissipation, magnetic 
reconnection or neutrino processes (capture on free nu- 
cleons, scattering on electrons and positrons or neutrino- 
antineutrino annihilation). Cooling occurs through neu- 
trino emission by nucleons and annihilation of electron- 
positron pairs. A detailed description of all these processes 
is beyond the scope of this paper and we have rather con- 
sidered two limiting cases in a very simplified way. 
When the disk is optically thin to neutrinos we adopt a 
uniform heating (per unit mass) qh along the field line and 
limit the cooling to neutrino emission by nucleons i.e. 



is veryi close to Vl ^except naturally for 6 = yir for which m * The cooling by the reyerse reactions ( neutrin0 em i s . 
00). The difference is typically less than 1% 6v6n for — 



sion by nucleons) is given by 



q e N 



= 1.4 10 i8 T b 



erg.g L .s 1 



(32) 



where T is the local temperature in the wind. Pauli block- 
ing effects for electrons have been neglected since 77 ~ 0.1 
everywhere in our wind solutions except in the vicinity 
of the disk. To compute the heating rate due to neutrino 
scattering on relativistic electrons and positrons we use 
the expression given by Herant et al. (1994 ) adapted to 
the disk geometry 

"4 



3.6 1 24 — f (T v -T)Tt 
P J 

disk 



dQ erg.g .s 



(33) 



The heating rate by neutrino-antincutrino annihilation is 
obtained from a double integral over the disk surface 



q = qh - QT 6 



(28) 



Although this is not s trictly satisfied in realistic m odels 
([Popham ct al. 1999t [VlacFadyen fc Woosley 1999j) , we 
require that q — in the plane of the disk which fixes qh 



disk disk 



(4) , (4.Y 



(1 — cosa)' 



-i (h I <M t„. I dn' i' Ve y (1-cosa) 



disk disk 



qh = QTq 



(29) 



(34) 



where Tp is the disk temperature at the foot of the line 
given by Eq. ( 8). 



where Qi and Q2 are two constants given in R.uffcrt 



et al. (1997 ) . The prime quantities correspond to a surface 
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Eq. (j28|). We obtain the mass loss rate rh in a classical 
way by integrating inward from the sonic point down to 
the disk surface. We fix trial values of T s and rj s at the 
sonic point from which we get v — v s from Eq. ( p5[ ) and 
the position y s from Eq. (p^). To start the integration 
we need the derivatives of v, T and 77 at y — y s . These 
derivatives cannot be directly calculated from Eq. (24) 
because F\ — F2 = F3 = A = at the sonic point. We 
instead use l'Hopital's rule which allows us to write three 
algebraic, second-order equations for the three derivatives 
and we solve these equations with a Newton-Raphson 
technique. Once v, T, rj and their derivatives have been 
determined at the sonic point the inward int egration 
can be started. In agreement with the results of Duncan 



Fig. 3. Wind solution for an optically thin disk with uni- 
form heating. The field line is anchored at r* = 4r/j and 
makes an angle 9 — 85° with the disk. The temperature 
at the foot of the line is Td — 2 MeV. In the velocity plot 
the dotted line shows the local sound velocity. The sonic 
point is located at y s — 2.174. The dashed line in the tem- 
perature plot corresponds to the approximate analytical 
solution (Eq. (47)). The electron degeneracy parameter 
r] = He/kT is close to 0.1 except at the vicinity of the 
disk. 

clement dS' whose neutrinos interact with those emitted 
by a surface element dS, a being the interaction angle. 
When the neutrino intensities and average energies are 
expressed as a function of the neutrinosphere temperature 
T v Eq. (Hf) becomes 

6 10 22 J dfl T^J dfl' T' V A {T V + T' u ){l-cosaf 




disk 



disk 



+ 6.7 io 20 | dn 1* J m' r,: 4 ^±^(i-cos«) 1 

disk disk / 

erg.g^.s -1 (35) 
Finally, the cooling rate from the annihilation of electron- 



positron pairs is given by (Herant et al. 1994) 



Qe+c 



rp9 

1.5 10 25 — er 



(36) 



4. The mass loss rate 

4.1. Numerical solution for an optically thin disk 

We first solve the flow equations for an optically thin disk 
and in the case of uniform heating, i.e. with q given by 



et al. (1986 ) for neutrino-driven winds in proto- neutron 
stars we observe that at some position y = y* the 
velocity begins to fall off rapidly while the temperature 
reaches a maximum T max < T d {r) (Td(r) being the disk 
temperature at radius r). We then adjust the values of T s 
and r) s with the requirement that should be as close as 
possible to and T max to T d (r). 

We first constructed a reference model where we follow 
the wind along a field line attached at r = r* =4rj. The 
disk temperature at r = r* is Td — 2 MeV. The line makes 
an angle 9 = 85° with the disk and the derivative 4^ = 0. 
The mass of the black hole is Mbh = 2.5 Mq. The results 
for this reference model are shown in Fig. |3[ The sonic 
point is located at y s = 2.174 slightly below y\ = 2.182 
where 7 = 0. The temperature and density at the sonic 
point are T s = 0.203 MeV and p s = 1.2 10 5 g.cm" 3 
which correspond to a sound velocity v s /c = 5.44 10" 2 
or v s = 16300 km.s" 1 and a mass loss rate rh = 2.3 10 14 
g.cm _2 .s _1 . The degeneracy parameter 77 remains practi- 
cally constant (77 ~ 0.1) from y = 0.5 to the sonic point, 
a property which will be used to construct the analytical 
solutions in Sect, [ll] below. Neutrino cooling is efficient 
close to the disk, up to y ^ 0.5. We have tested the effect 
of a different field line geometry with 9' = W- ^ 0. We 
considered two cases, 9' = — 2° /rh and — 5°/rft, for which 
we respectively obtain m = 2.5 and 2.8 10 14 g.cm" 2 .s _1 . 

We then varied the mass of the black hole Mbh, the 
disk temperature Td, the inclination angle 9 and the posi- 
tion r of the foot of the line to see how these parameters 
affect the mass loss rate. The results are shown in Fig. ^ 
where we have plotted rh when Mbh, Td, 9 and x — r/r^ 
are varied separately while the other three quantities are 
maintained at fixed values (choosen to be those of the re- 
ference model: M B h = 2.5 M , T d = 2 MeV, 9 = 85° 
and x = 4). The mass loss rate appears to be nearly pro- 
portional to the mass of the black hole (Fig. ||a). The 
dependence of m on the temperature Td at the foot of the 
field line is much more spectacular since we get rh ex Tj° 
(Fig. |Jb) in agreement with the results for neutrino-driven 
winds in neutron stars ( Duncan et al. 1986| ). The mass loss 
rate also sharply increases when the angle between the 
field line and the disk is reduced (Fig. ||c) . Below 9 ~ 77° 
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it becomes more and more difficult to construct wind so- 
lutions with the required accuracy. This can be related to 
the value of the Bcrnouilli function 



B = -v 2 + h + $ - $1 



(37) 



(h being the specific enthalpy and $i the potential at y = 
?/i ) which is positive in the plane of the disk for 9 < 77°. 
The initial thermal energy is then sufficient to allow the 
gas to escape even in the absence of additional heating. 
Finally, the mass loss rate approximately increases as r 
(Fig. [|d) when the potential well becomes shallower at 
larger radial distances. 

As long as the Bernouilli function is not too close to zero 
the mass, temperature and geometrical dependence can 
be separated in m to yield the general expression 



m(x) « 2.3 10 14 a<bh 



T d (x) 



2 MeV 



10 



f[x,8(x)} g.cm 2 .s 1 



(38) 



where / is a geometrical function which satisfies 
/(4,85°) = f. The mass loss rate has been represented 
in Fig. U as a function of x for /ibh = 1, T^{x) given by 
Eq. (||) and different field geometries: constant inclination 
angles 9 = 80, 85 and 89° or 9 linearly decreasing from 90° 
at x = 3, i.e. 9 = 90° - \{x - 3) (with A = 0.5, 1.0, 1.5 or 
2°). At constant 9 the mass loss rate varies rather slowly 
with x while for A = 1.5 or 2° it rapidly increases at large 
radial distances where 9 < 80°. 



4.2. Numerical solution for an optically thick disk 

If the disk is optically thick the heating and cooling 
sources due to neutrinos are specified by Eq. (|3l|-|3^). 
Again, wind solutions are found by inward integration 
from the sonic point down to the disk. A reference model 
is constructed with the same black hole mass and field 
geometry (A/bh = 2.5 Mq and field line attached at 
r = r* = 4 r/j and making an angle 9 = 85° with the 
disk). The disk is supposed to be optically thick to neu- 
trinos from nn = 3 rh to r out = 10 r/,. The temperature 
of the neutrinosphere in this reference model is T v = 2 
MeV and does not vary with radius which allows a sim- 
ple calculation of the geometric integrals appearing in 
the neutrino heating terms. The results for this referen- 
ce model are shown in Fig. |6| The sonic point is located 
at y s = 2.175. The temperature and density at the sonic 
point are T s = 0.132 MeV and p s = 3.23 10 4 g.cm" 3 
which correspond to a sound velocity v s /c — 4.45 10~ 2 or 




Fig. 5. Mass loss rate per unit surface as a function of 
x = r/fh for an optically thin disk with uniform heating. 
Upper panel: field lines with constant inclination angles 
9 = 80, 85 and 89°; lower panel: field lines with decreasing 
inclination angle 9 — 90° — X(x — 3) with A = 0.5, 1, 1.5 
and 2° (curves labelled (a) to (d)). 




13350 km.s 1 and a mass loss rate m = 5.1 10 13 Fig. 6. Same as Fig. ||for an optically thick disk. 



g.cm~ 2 .s -1 . The neutrino heating and cooling terms are 
detailed in Fig. [?]. The major contribution to the heating 
comes from neutrino captures on nucleons while neutrino 
emission by nucleons and annihilation of electron-positron 
pairs have comparable effects on the cooling. We have then 
abandoned the assumption of constant neutrino tempera- 
ture which was essentially made for a rapid calculation of 



the integrals in the neutrino heating terms. When T u de- 
pends on the radial distance in the disk, the calculation 
of q V y becomes very time consuming. Since <jw is not the 
dominant neutrino heating term as long as T v <> 10 MeV 
we have neglected its contribution for the remainder of 
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2.5 5 7.5 10 1 2 

M BH (M e ) T d (MeV) 




75 80 85 90 4 6 8 10 



00 x = r / r h 

Fig. 4. Mass loss rate per unit surface for an optically thin disk with uniform heating as a function of (a) black hole 
mass Mbh, (b) disk temperature at the foot of the line Td, (c) inclination angle 9 and (d) radial distance to the black 
hole. When one quantity is varied the three others are maintained at their reference values: Mbh = 2.5 Mq, Td = 2 
MeV, 9 = 85° and r — Art- The dotted lines in (a), (b) and (d) have respective slopes 1, 10 and 4.4. 



this paper. For example, in our reference model the mass 
loss rate is decreased from 5.1 to 4.9 10 13 g.cm _2 .s _1 if 
q v p is not included (i.e. a reduction of 4.5%). Without q V y, 
a non constant neutrino temperature can be easily imple- 
mented. We adopted for T„(r) relation (|5|) with T* = 2 
MeV and the resulting mass loss rate is then m = 3.8 10 13 

— 2 — 1 

g.cm .s . 

When Mbh, T v , 9 or the position r of the foot of the line 
are varied, the mass loss rate behaves as in the optically 
thin case 



m(x) « 3.8 10 /xbh 



T v {x) 



1 io 



2 MeV 



f[x,9(x)] 



g.cm 2 .s 1 



3.8 10 13 /i BH 



2 MeV 



(t) 



15/2 / i _ ./m 



10 



5/2 



1 - 



(39) 



The mass loss rate has been represented in Fig. |8| for T* = 
2 MeV, r* = 4 rft, ri„ =3^, r ou t = 10 and the same 



field geometries already considered in Sect. [11 
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0) 



10 1 
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VV 



vN 



ve 



Fig. 7. Neutrino heating and cooling contributions. The 
long and short dashed lines respectively represent indi- 
vidual heating and cooling processes; vN: neutrino cap- 
tures on nucleons; ve: neutrino captures on electrons; vv: 
neutrino-antineutrino annihilation; eN: neutrino emission 
by nucleons; e + e~: annihilation of electron-positron pairs. 
The full thin lines represent the total of the heating (vN 
+ ve + vv) and cooling (eN + e + e~) processes. The full 
thick line is the sum of all contributions. 



4.3. Analytical solution 

To obtain an analytical expression for the mass loss rate 
we have simplified the original wind problem by making 
several additional assumptions. We have first considered 
that wind material is not too far from hydrostatic equili- 
brium even at the sonic point (in practice, the ratio 
becomes larger than unity at v "~ v ~ 5%). We have also 
supposed that the pressure is dominated by the contribu- 
tion of relativistic particles, which is equivalent to 77 <C 1 
since 



Pn 

Pr 



0.5 



Y, 



(40) 



In practice 77 ~ 0.1 in our numerical solutions, except near 
the disk surface but we have nevertheless adopted P ~ 
P r oc T 4 everywhere (we take P r = dT 4 with 5 being 
given by Eq. ([IJ) with 77 = 0.1). Finally, we have located 
the sonic point at y s — yi where j(y) — 0. We then write 
the temperature 



T = T S 



1/3 



(1 + r) 



(41) 



where the r function is zero at the sonic point. From the 
sonic point down to the disk, as i] cx p/T 3 first remains 
nearly constant (see Fig. || and ^) r stays close to zero. 



E 
u 

\ 

cn 

•E 



I0 16 
I0 15 
I0 14 

I0 12 



n 1 1 1 r 



~i 1 1 1 1 r 



H — h 



— I — I — I — ^ 


— I — I — I — 








(d! 




Fig. 8. Same as Fig. |]for an optically thick disk. 

It then becomes negative when r\ increases in the vicinity 
of the disk. With these assuptions the wind equations can 
be rewritten 



pv s = m 

P dlogT d$ 

4 = 7 r = — 

p dy dy 

d\og (1 + t) qr 



dy 



3v 



ps qr 
3m 



The solution of Eq. (43) gives 



Tn-T = 



1 



1 



4^9 J (1 + r) 3 dy 



dy 



(42) 
(43) 

(44) 
(45) 



where 9 = T s 3 / p s and To is the temperature in the plane 
of the disk. Taking r = everywhere leads to the appro- 
ximate solution 



T -T 



1 



4&e 



My) - *o] 



(46) 



with $0 = Q(y = 0). If we moreover suppose that T s To 
we get 



T $ s -$(y) 

t(y) = — ~ — ^ 

Ky > To $ s - $ 



and 



6 



AaTo ~ 4aTo 



(47) 



(48) 



with <j) s — $(j/ s ) = &(yi)- I n Fig- p| th e approximate value 
of the temperature given by Eq. (|47| ) is compared to the 
exact solution obtained in Sect. 4.1. It can be seen that 



the agreement is quite satisfactory for y < 1. We then 
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transform Eq. ( [44| ) using relation (f4l|), the definition (| 
of t{y) and writing q{y) as 



q (y) = q n [g(v) - t 6 (v)} 



(49) 



where g(y) = 1 for an optically thin disk with uniform 
heating. In an optically thick disk we only consider the 
heating due to neutrino captures on nucleons and g(y) is 
the geometrical integral appearing in Eq. (|3l|). We obtain 

2 02 d(l+r) 6 _ s{y)Qr rrB r+2 



dy 



3m 



T${t 2 (y) [g(y)-t 6 (y)}} (50) 



Integrating this equation from the disk to the sonic point 
yields 



2ae 2 [l-(l + r ) 61 " 
where the integral 



A$ 2 -/i 2 



s(y)t*(y)\g(y)-t 6 (y)] dy 



can be directly computed from Eq. p^) , ([ 
obtain Eq. (|l|) we have used the value 



45T 4 



(1 + t \ a = — — = !) = -^L 
' p Q p A$ A$ 

all quantities with a zero subscript being taken in the 
plane of the disk. The analytical formula (|Tj) reproduces 
the Tq° dependence of to and diverges when ho — A$, 
i.e. when the Bernouilli function is equal to zero at the 
disk surface. The mass loss rate is also proportional to 
the mass of the black hole at constant x — r/r^. We have 
tested the accuracy of the analytical expression by com- 
puting to for the reference model with uniform heating 
(optically thin disk). We get to = 3.9 10 14 g.cm~ 2 .s~ 4 in 
reasonable agreement with the numerical result obtained 



(52) 

2) and ©. To 
P of O and 

(53) 



in Sect. 4.1 



5. Estimates of the wind Lorentz factor 

In the simple model considered in this paper we do not fol- 
low the acceleration of the wind beyond the sonic point up 
to relativistic velocities. Therefore, the terminal Lorentz 
factor cannot be obtained in a self-consistent way. We sim- 
ply expect that an average value of the terminal Lorentz 
factor will be given by 

L - (54) 



r = 



Mc 2 

where L w is the power injected into the wind and M is 
the total mass loss from the disk 



M 



to lixrdr 



(55) 



The mass loss per unit surface of the disk has been ob- 
tained in the previous section 



TO(x) = TO £tBH 



T d (x) 



2 MeV 



f[x,6(x)] 



(56) 



so that 



M = 3.4 10 12 too a*| h 





' T d (x) ' 




2 MeV 



1 10 



f[x,6{x)]xdx (57) 



We have computed M for x- ln = 3 and x out = 10, both 
for optically thin and optically thick disks. In optically 
thin disks the temperature is given by Eq. (^) and too = 
2.33 10 14 g.cm" 2 ^- 1 . We have then 



M = 5.1 10 28 /ibh 
with 



0.01 



gco 



(51) F„ 



f[x,6(x)] 



dx 



(58) 



(59) 



For an optically thick disk the temperature of the neutri- 
nosphere is given by Eq. (^) and too = 3.8 10 13 g.cm _2 .s _1 
so that 



M = 7.9 10 2 VIh f 1 



2 MeV J 



g.s 



(60) 



with 



F„ 



X 



0.5 



5/2 



f[x,6{x))dx (61) 



We now introduce the ratio (3 of the wind power to the 
disk neutrino luminosity 



L W = (3L V = (3J -<rT*(r) 2nrdr 
which yields the following estimate of T* 



F* — 2 /igj/ 



L$2 

IT 



1/4 



MeV 



(62) 



(63) 



where L 52 = L w /10 52 erg.s -1 and /3g = (3/5. A large (3 
corresponds to a disk able to transfer most of its energy 
in the Poynting flux with limited dissipation and heating. 
From Eq. (p8|) and (63) we finally get 



M = 7.9 10-Mbh(^ 



-52 



5/2 



gco 



g.S 



(64) 



We are now in a position to estimate the average 
Lorentz factor of the wind using Eq. (|54|). The results 
in the optically thin and optically thick cases are given 
respectively by 



T = -^^11^52 ( o^x ) 2 opt. thin case 
= -^-/^bh (3^ 2 opt. thick case 



(65) 



and the values of the geometric integrals are listed in 
Table |. 



12 



F. Daigne & R. Mochkovitch: Baryonic pollution in GRBs: the case of a magnetically driven wind 



Geometry 


Opt. thin disk 


Opt. thick disk 




fgco 


r 


-TgCO 


r 


9 = 80° 


9.0 


24 


1.1 


130 


= 85° 


1.1 


200 


9.2 lO -2 


1500 


= 89° 


0.36 


610 


2.6 HT 2 


5400 


<9 = 90°-0.5°(x-3) 


0.55 


400 


4.0 10-* 


3500 


= 9O°-l.°(a;-3) 


1.3 


170 


9.1 10" 2 


1500 


= 90° -1.5° (a? -3) 


> 5.8 


< 38 


0.42 


340 


<9 = 90°-2.°(:r-3) 


> 200 


« 1 


> 3.3 


< 43 



Table 1. Values of the geometric integral F goo for an opti- 
cally thin and an optically thick disk and the field geome- 
tries considered in Sect. 4. The corresponding values of f 
are computed using Eq. |65| with /ibh = 1, £52 = 1 and 
a = 0.01 in the optically thin case and /ibh = 1, £52 = 1 
and j3 = 5 in the optically thick case. 



6. Discussion and conclusions 

From Eq. ( |65| ) and the values of the geometric integrals 
m Table | it appears that the wind emitted from the 
inner disk can reach large Lorentz factors (r > 100) only 
under quite restrictive conditions on the disk temperature 
and field geometry. The disk temperature depends on the 
viscosity parameter a (in the optically thin case) and on 
the value of (3 (in the optically thick case). These two 
quantities which are quite uncertain unfortunately enter 
the expression of T with respective powers 2 and 2.5. 
The extreme sensitivity of T to the disk temperature also 
amplifies all errors and uncertainties in the evaluation of 
Trf. For example, the simple analytical expression of Td 
(Eq. 8) used in the optically thin case does not depend on 
the accretion rate to the black hole while detailed calcu- 



lations sho w that slightly increases with M (Popham 
et al. 1999| ). A consequence of Eq. (8) is that f oc L but 



just the opposite behavior (decreasing T with increasing 
L) would be expected if Td oc M x with x > 0.1. It is 
moreover even not clear that the a-prescription remains 
appropriate in the context of strongly magnetized disks. 

The Lorentz factor is also extremely sensitive to the 
field geometry: to escape from the disk, the material has 
to be heated so that its energy becomes large enough 
to cross the potential well. The height of this barrier 
decreases rapidly when the field lines are more inclined, 
leading to an increasing mass flux. As shown by our 
analytical solution of the wind equations, there is even 
a critical inclination where the mass flux diverges. This 
happens when the potential barrier is so shallow that 
the initial enthalpy of the material in the disk allows 
it to freely escape along the field lines, even without 
additional heating. For this reason quasi vertical field 
lines are required to prevent baryonic pollution from 
growing dramatically and slight changes of the inclination 
angle lead to large variations of the Lorentz factor. 
Between 9 = 90 and 80°, f approximately behaves as 
8 3Q ! A reduction by 1° of the inclination angle therefore 



decreases T by a factor of about 1.5. 

Our incomplete description of the black hole and its 
environment prevents us from making any accurate pre- 
diction of the wind Lorentz factor. The numerical values 
of the mass loss rate given above have been obtained with 
a large set of simplifying assumptions. We present them to 
illustrate general tendencies which we believe are robust 
and can be used to evaluate the ability of realistic models 
to produce relativistic outflows. When the disk is thick 
or slim the temperature distribution and the heating or 
cooling processes will be different from the assumptions 
we have made. Similarly, the disk + black hole magneto- 
sphere can be expected to have a very complex geometry 
and to be rapidly variable. The results of our toy model - 
extreme sensitivity of the mass loss rate to the disk tem- 
perature (rh oc Tj°) and field geometry - indicate that 
high terminal Lorentz factors can be reached only under 
severe constraints: disk temperature Td not largely exceed- 
ing 2 MeV and presence of at least a few field lines pointing 
directly away from the disk in the vertical direction. 

If such conditions can be satisfied so that the wind 
can indeed become relativistic, it is then easy to under- 
stand that its baryonic load and hence its Lorentz fac- 
tor can strongly vary on short time scales as a result 
of fluctuations of the disk temperature or field geometry. 
Conversely, if the outflow remains non relativistic (r ~ 1) 
because the disk is too hot or the field lines deviate too 
much from the vertical, the burst must be produced by the 
Blandford-Znajek effect alone with no contribution from 
accretion energy. The dense wind emitted from the disk 
can then have both a beneficial and negative effect on the 
central relativistic jet. It can probably help to confine the 
jet but can also represent a risk of baryonic pollution via 
Kelvin- Helmholtz instabilities or magnetic reconnection at 
the jet-wind interface. 

Acknowledgements. We thank the anonymous referee for care- 
ful reading of the manuscript and thoughtful comments. 
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